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Abstract. We consider new concepts of entropy and pressure for stationary 
systems acting on density matrices which generalize the usual ones in Ergodic 
Theory. Part of our work is to justify why the definitions and results we 
describe here are natural generalizations of the classical concepts of Thermo- 
dynamic Formalism (in the sense of R. Bowen, Y. Sinai and D. Ruelle). It is 
well-known that the concept of density operator should replace the concept of 
measure for the cases in which we consider a quantum formalism. 

We consider the operator A acting on the space of density matrices Mn 
over a finite TV-dimensional complex Hilbert space 



where Wi and Vi, i = 1, 2, , k are linear operators in this Hilbert space. In 

some sense this operator is a version of an Iterated Function System (IFS). 
Namely, the Vi (.) V* =: Fj(.), i = 1,2, ...,k, play the role of the inverse 
branches (i.e., the dynamics on the configuration space of density matrices) 
and the Wi play the role of the weights one can consider on the IFS. In this 
way a family W := {Wi}i—\ ... fc determines a Quantum Iterated Function 
System (QIFS). 

We also present some estimates related to the Holevo bound. 
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In this work we investigate a generalization of the classical Thermodynamic 
Formalism (in the sense of Bowen, Sinai and Ruelle) for the setting of density 
matrices. We consider the operator A acting on the space of density matrices M. n 
over a finite iV-dimensional complex Hilbert space 



where Wi and Vi, i = 1,2, ... ,k are linear operators in this Hilbert space. Note 
that A is not a linear operator. This operator can be seen as a version of an Iterated 
Function System (IFS). Namely, the V % (.) V* =: F 4 (.), i = 1, 2, . . . , k, play the role 
of the inverse branches (i.e., the dynamics on the configuration space of density 
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matrices p) and the W% play the role of the weights one can consider on the IFS. 
We suppose that for all p we have that 2<=i t r (WipW*) = 1. Note that such trace 
preserving condition, for any normalized operator p (that is, with tr(p) = 1), is 
equivalent to the explicit condition J2i W*Wi = I. We say that A is a normalized 
operator. 

A family W :— {Wi}j=i,...,fc determines a Quantum Iterated Function System 
(QIFS) T w , 

T w = {M N ,F l ,W i } i = h ..., k 

Basic references on QIFS are and [TB]- We want to consider a new concept 
of entropy for stationary systems acting on density matrices which generalizes the 
usual one in Ergodic Theory. In our setting the Vi, i — 1,2, ... ,k are fixed (i.e. 
the dynamics of the inverse branches is fixed in the beginning) and we consider 
the different families Wi, i = 1,2, ...,k, (also with the attached corresponding 
eigendensity matrix pw) as possible Jacobians of stationary probabilities. 

Given a normalized family Wi, i — 1,2, ... ,k, a natural definition of entropy is 
given by 
(2) 

hv{w) = " § Hv Pw v;) § »{W"V* W J ) ^ ( tr(V iPw V*) , 

where pw denotes the barycenter of the unique invariant, attractive measure for the 
Markov operator V associated to Tw ■ We show that this generalizes the entropy 
of a Markov System. 

We also want to present a concept of pressure for stationary systems acting on 
density matrices which generalizes the usual one in Ergodic Theory. In addition 
to the dynamics obtained by the Vi, which are fixed, a family of potentials Hi, 
i = 1, 2, . . . h induces a kind of Ruelle operator given by 

A' 

(3) C H (p):=Y,tr(H iP H*)VipV* 

i=l 

We show that such operator admits an eigenvalue j3 and an associated eigenstate 
P/3, that is, one satisfying Cn{pp) = Pp- 

The natural generalization of the concept of pressure for a family Hi, i = 
1,2, ... ,k is the problem of maximizing, on the possible normalized families Wi, 
i = 1, 2, . . . k, the expression 

fe 

(4) h v (W) + log (tr{H 3P[j H*)tr{V jP pV*))tr{W p w W*) 

3=1 

We show a relation between the eigendensity matrix pp for the Ruelle operator and 
the set of Wi, i = 1, 2, ... k, which maximizes pressure. In the particular case that 
each of the Vi is unitary, i = 1, 2, ... k, the maximum value is log /3. 

Our work is inspired by the results presented in pi)] and [21] . We would like to 
thank these authors for supplying us with the corresponding references. 
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It is well-known that completely positive mappings (operators) acting on density 
matrices are of great importance in Quantum Computing. These operators can be 
written in the Stinespring-Kraus form (see section [T2"]) . Also a nice exposition on 
the interplay of Ergodic Theory and Quantum Information is presented in [4] . 

The initial part of our work aims to present some of the definitions and con- 
cepts that are not very well-known (at least for the general audience of people in 
Dynamical Systems), in a systematic way. We present the main basic definitions 
which are necessary to understand the theory. However, we do not have the inten- 
tion of exhausting what is already known. We believe that the theoretical results 
presented here can be useful as a general tool to understand problems in Quantum 
Computing. 

Several examples are presented in the text. We believe that this will help the 
reader to understand some of the main issues of the theory. In order to simplify 
the notation we will present most of our results for the case of matrices of order 2. 

In sections [2] and[3] we present some basic definitions, examples and we show some 
preliminary relations of our setting to the classical Thermodynamic Formalism. In 
section [4] we present an eigenvalue problem for non-normalized Ruelle operators 
which will be required later. Some properties and concepts about density matrices 
and Ruelle operators are presented in sections [6] and [7] Sections [8] and [9] are 
dedicated to the introduction of some different kinds of entropy that were already 
known but do not have a stationary character. In section [10] we introduce the 
concept of stationary entropy for measures defined on the set of density matrices. In 
section[TT]we compare this definition with the usual one for Markov Chains. Section 
IT2l is dedicated to motivate the interest on pressure and the capacity-cost function. 
Section [131 [14] [15] and [16] are dedicated to the presentation of our main results on 
pressure, important inequalities, examples and its relation with the classical theory 
of Thermodynamic Formalism. 

In [T] we present a general exposition (describing the setting we consider here) 
where we omit proofs, but provide many examples. We believe that paper will help 
to complement the present paper for the reader which is a newcomer in the area. 
We also present there some basic results concerning the discrete Wigner measure. 

In [5] we propose a different concept of entropy which is also a generalization 
of the classical one. We also describe some properties of the Quantum Stochastic 
Process associated to the Quantum Iterated Function System. 

This work is part of the thesis dissertation of C. F. Lardizabal in Prog. Pos-Grad. 
Mat. UFRGS (Brazil). 

2. Basic definitions 

Let Mat(C) the set of complex matrices of order n. If p 6 Mjv(C) then p* denotes 
the transpose conjugate of p. A state (or vector) in C™ will be denoted by -0 or 
and the associated projection will be written \ijj)(t/)\. Define 

U N := {p G M N (C) : p* = p} 
VH N := {peH N : > 0,W> e C N } 
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M N ■= {peVHN ■ tr(p) = 1} 

V N := {p e H N : p = e C", = 1}, 

the space of hermitian, positive, density operators and pure states, respectively. 
Density operators are also called mixed states. If a quantum system can be in one 
of the states {ipi, . . . , ipk} then a mixed state p will be written as 

k 

(5) p = ^2pi\i>i){i>i\ 

i=l 

where the pi are positive numbers with J2iP* = 1- 



Definition 1. Let Fi : Mn — > Mn , Pi ■ Mn — > [0,1], i = 1, • • ■ , k and such that 
Y,iPi(p) =1- We call 

(6) F N = {Mn,Fi, pi : i = 1, ...,k} 

a Quantum Iterated Function System (QIFS). 

Definition 2. A QIFS is homogeneous if Pi and FiPi are affine mappings, i — 
l,...,fc. 



Suppose that the QIFS considered is such that there are V% and W, linear maps, 
i=l,...,k, with J2i=i W*Wt = I such that 

(?) FliP) trJV~p~V*) 

and 

(8) p i (p) = tr(W iP W*) 

Then we have that a QIFS is homogeneous if Vi=Wi, i = 1, . . . , k. Now we can 
define a Markov operator V : M 1 (Mn) -> A1 1 (A^jv), 

(9) (V M )0B)=W Pi{p)dp{p), 

4 =i^r 1 (s) 

where A^ 1 (A^at) denotes the space of probability measures over Mn- We also 
define A : Mn ~^ Mn, 

k 

(10) A(p) := 

i=l 

The operator defined above has no counterpart in the classical Thermodynamic 
Formalism. We will also consider the operator defined on the space of density 
matrices p, 

k 

(11) C(p) = J2<n(p)VipV*. 

i—k 
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If for all p we have J2i=k 9i(p) = 1, we say the operator is normalized. We are also 
interested in the non-normalized case. If the QIFS is homogeneous, then 

(12) k(p)=J2VipV* 

i 

Theorem 1. |21j A mixed state p is A-invariant if and only if 

(13) p= pdp{p), 

J Mn 

for some V -invariant measure p. 

We recall the definition of the integral above in section 

In order to define hyperbolic QIFS, one has to define a distance on the space of 
mixed states. For instance, we could choose one of the following: 

(14) D 1 (p 1 ,p 2 ) = y/trftpj. - p 2 f] 

(15) D 2 ( Pl ,p 2 ) = try/( Pl -p 2 ) 2 

(16) D 3 (p uP2 ) = VV-M^W 2 ) 1 / 2 ]}, 

the Hilbert-Schmidt, trace, and Bures distances, respectively. Such metrics generate 
the same topology on A4n- Considering the space of mixed states with one of those 
metrics we can use a definition of hyperbolicity similar to the one used for IFS. That 
is, we say a QIFS is hyperbolic if the quantum maps F% are contractions with 
respect to one of the distances on jMw and if the maps pi are Holder-continuous 
and positive, see for instance [16] . 



Proposition 1. If a QIFS {IP is homogeneous and hyperbolic then the associated 
Markov operator admits a unique invariant measure p. Such invariant measure 
determines a unique A-invariant state p £ Mn, given by (TO)). 

See [16], [21] for the proof. 

3. Examples of QIFS 

Example 1. ft = M N , k = 2, Pl = p 2 = 1/2, G x {p) = U lP U?, G 2 {p) = U 2 pU^. 
The normalized identity matrix = I/N is A-invariant, for any choice of unitary 
Ui and U 2 . Note that we can write 

(17) p* = / pdp(p) 

where the measure p, uniformly distributed over Vn (the Fubini-Study metric), is 
V-invariant. 
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We recall that a mapping A is completely positive (CP) if A ® I is positive 
for any extension of the Hilbert space considered Hn <S>He- We know that 

every CP mapping which is trace-preserving can be represented (in a nonvmiquc 
way) in the Stinespring-Kraus form 

k k 

(is) K( P ) = Y J v 3 pv;, £v7^- = i, 

j=l ] = 1 

where the Vi are linear operators. Moreover if we have ^2j = iVjV* — I, then 
A(I/N) = I/N. This is the case if each of the V\ are normal. 

We call a unitary trace-preserving CP map a bistochastic map. An example 
of such a mapping is 

k 

(19) Au(p) = Y,PiUipUi, 

z=l 

where the Ui are unitary operators and YliPi = 1- Note that if we write i*i(p) = 
UipU* , then example[T]is part of this class of operators. For such operators we have 
that /J* is an invariant state for Kjj and also that S Pt is invariant for the Markov 
operator Py induced by this QIFS. 

We will present a simple example of the kind of problems we are interested 
here, namely eigenvalues and eigendensity matrices. Let Wjv be a Hilbert space 
of dimension N. As before, let A4n be the space of density operators on Hn- A 
natural problem is to find fixed points for A : M n ~ > -M n > 

k 

(20) A(p) = J2 V iP V * 

In order to simplify our notation we fix N = 2 and k = 2. Let 

v = ( Vl V2 V v 2 = ( Wl W2 V P =( 1 p > 

\ v 3 v A J \ w 3 Wi J \ p 2 Pa 

where V\ and V% are invertible and p is a density operator. We would like to find 
p such that 

(21) V lP V* + V 2P V* = p. 
Example 2. Let 



where fc, I e M, p £ (0, 1). Then V^Vi + V 2 *V 2 = I. A simple calculation shows that 
p 2 = 0, and then 

q 
1-q 

is invariant to A(p) = V\pV( + V 2 pV 2 * , for q G (0, 1). 





P = 
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Now we make a few considerations about the Ruelle operator L denned before. 
In particular, we show that Perron's classic eigenvalue problem is a particular case 
of the problem for the operator C acting on matrices. Let 

r. 



f Poo 




, v 2 = ( 


f 


Vox \ 


V o 


J 


. o 


) 



V3={ o o\ n _fo X „_( P1 p 2 



PlO / V pxi J \ P3 Pi 

Define 

4 

(22) C(p) = J2<n(pWipV* 

i=l 

We have that C(p) — p implies p 2 = and 

(23) api + bp4 = pi 

(24) epi + dpi = pi 
where 

a = qiplo, b = q 2 p : 1 , c = q 3 p 2 10 , d = q 4 pl x 
Solving ([23)1 and in terms of p\ gives 



b l-d 

Pi = -, Pa, Pi = Pi 

1 — a c 



that is, 
(25) 



1 -a 



which is a restriction over the qi. For simplicity we assume here that the qi are 

constant. One can show that 

(26) 

P = 
Now let 



92Poi-<?lPo +l I _ I l-q4Pii+93Pio 

o *zgy$ J I o , — T i 2 

92Poi-9lPo +l / \ 1-94PI 1 +93PI 



be a column-stochastic matrix. Let ir = (771,772) such that P77 = 77. Then 



Poo P01 
P10 P11 



(27) rr=( ™ , ) 

P01 - Poo + 1 poi - Poo + 1 

Comparing (|27|) and (|26|) suggests that we should fix 

1 1 1 1 

(28) gi = , q 2 = , q 3 = , qi = 

Poo P01 P10 P11 
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Then the nonzero entries of p are equal to the entries of 7r and therefore we associate 
the fixed point of P to the fixed point of some £ in a natural way. But note that 
such a choice of qi is not unique, because 

(oa\ n 1 Z qip m „ 1 - gsPioPoi 

l^yj Q2 = , <?4 = 2 ' 

PoiPio Pu 
for any qi , q 3 also produces p with nonzero coordinates equal to the coordinates of 
7r. We also note that the above calculations can be made by taking the Vi matrices 
with nonzero entries equal to y/pij instead of Pij . 

Now we consider the following problem. Let 

f h 00 \ f h 01 \ _ f 

Vl " I o o ) ' V2 { o ) ' v * Uo o 

Define 



/in / ^ \ p 3 P4 



(so) = 

i=l 

where % e R. Assume that fey € R, so we want to obtain A such that C{p) = Xp, 
A 7^ 0, and A is the largest eigenvalue. With a few calculations we obtain p 2 = pz = 
0, 

lihloPi +q2hl 1 p i = Api 
<lz h wPi + 94^?iP4 = Ap 4 

that is, 

(31) api + bp 4 = Api 

(32) cpi + dp 4 = Ap 4 , 
with 



Therefore 



and 



a = <7ifyro; b = Q2hl 1 , c = q 3 hl , d = g 4 /i?i 



^ OA/ ^p 4 



/04 / y pi 
X-d b 



c X — a 
Solving for A, we obtain the eigenvalues 



where 



a + d (_a + d ^ \d - a) 2 + Abe 
A — dz — — dz 

2 2 2 2 
= \ (<?i^oo + <H h ii ± \J {q~J^ 1 - qihl ) 2 + 4q^ji%A) 1 

C= ^{d- a) 2 +4bc = xj (q A h 2 n - qih^) 2 + ^©fe^io 
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and the associated eigenfunctions 
P = 

But p\ + pi = 1 so we obtain 



0\ ^p 4 

Pi ) \ pi 



a-d±C 



n _ I a~d±C+2c 

p -\ 



a-d±C+2c 



(33) 
that is, 



qih (m -q4h 11 ±( 



-26 



(34) 







-2g 2 / I g 1 
9i' l oo- 2, ?2'ioi-'?4/ii 1 =FC 





9i /l oo _ 94' l ii ± C+ 2 93'iio 





a— d=pC 



a-2b-d=pC 





q 1 hl -2q 2 hl 1 -q i h , { 1 TQ 



Therefore we obtained that pi, p&, qi, , . . , X are implicit solutions for the set of 
equations (| 3 1 1) - (f32 |) . Recall that in this case we obtained pi = pz = 0. 

Now we consider the problem of finding the eigenvector associated to the domi- 
nant eigenvalue of H . The eigenvalues are 

1 



A = - [h m + h n ± \/(h 00 - /in) 2 +4h 01 h 10 
Then we can find v such that Hv = Xv from the set of equations 

(35) h QO vi + h iv 2 = Xvi 

(36) h w vi + huV2 = Xv 2 
which determine V\ , v 2 , A implicitly. Note that if we set 



(vr\ 1111 

(37) qi = , q 2 = , 93 = , 94 = 

Poo Pol PlO Pll 



we have that the set of equations (|31l) - (|32j) and (|35l) - (|36|) are the same. Hence we 
conclude that Perron's classic eigenvalue problem is a particular case of the problem 
for C acting on matrices. 



A different analysis in the quantum setting which is related to Perron's theorem 
is presented in [6]. 
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4. A THEOREM ON EIGENVALUES FOR THE RUELLE OPERATOR 

The following proposition is inspired in |18j . We say that a hermitian operator 
P : V — > V on a Hilbert space (V, (•)) is positive if (Pv,v) > 0, for all v E V, 
denoted P > 0. Consider the positive operator Cw,v ■ PT~Ln — > VHn-, 

k 

(38) Cw,v{p) ■■= Y.HWipWftVipVi*. 

i=l 

We point out that this operator is completely general. In an analogy with the 
classical case we can say it corresponds to the general Perron Theorem for positive 
matrices (having positive eigenvalues which can be bigger or smaller than one), by 
the other hand the setting described in [IB], [2J "basically" considers the analogous 
case of the Perron Theorem for stochastic matrices. 

We need a result in this form in order to better understand the Pressure problem 
which will be described later. 

Proposition 2. There exists p € Mjv and (3 > such that £w,v(p) = ftp. The 
value P is obtained explicitly: (3 = tr{Lw.y{p))- 

Proof Define C n : M n -> M n , 

r I \ £w,v{p+ i) 

The operator above is well defined. In fact, note that Cw,v(p), WjW* , VjV* are 
positive for all j. Then 



i i 

= Y J tr{W i pW* + -WiW*)tr{V iP V* + -V.V*) > 
* — ' n n 

> J2 tr (WipW*)tr{V iP V*) - tr(C W y) 

i 

We know that for any positive operator P ^ 0, if {vi, . . . ,«v} is a orthonormal 
base for Hn, then 

N 

tr(P)=Y,(Pvi,Vi) >0 

i=l 

Therefore, tr(Cw,v(p + r)) > 0, n > 1. Hence C n {p) is well defined. 

We know that Mn is compact and convex, so we can apply Schauder's theorem 
for each of the mappings £„, n > 1 and get p n € Mn such that 

Cn(Pn) = Pn => Cw.vipn H ) = PnPn, U > 1 

n 

where 

fi n := tr{C Wt v{p n + -)) 
n 
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By the compacity of Mn, we can choose a point p £ Mn which is limit of the 
sequence and then, by continuity, Cw,v(p) = ftp, where ft = tr{Lw,v{p))- 

Also, note that ft > 0, because if {vx, . . . , ujv} is a orthonormal base of Hn, 

N 

tr(£ wy (p)) = y^ y {Cw,v(p)Vi,Vi) > 0, 

i=l 

since Cw,v(p) is positive, and the inequality will be equal to zero if and only if 
£w,v(p) is the zero operator. Hence, we proved that there exists p e Mj/ and 
ft > such that Cw,v{p) — ftp- 

□ 

5. Vector integrals and barycenters 

We recall here a few basic definitions. For more details, see [16] and [21] , Let 
X be a metric space. Let (V, +, •) be a real vector space, and r a topology on V. 
We say that (V, +, •; r) is a topological vector space if it is Hausdorff and if the 
operations + and • are continuous. For instance, in the context of density matrices, 
we will consider V as the space of hermitian operators 'Hjv and X will be the space 
of density matrices Mn- 

Definition 3. Let [X, E) be a measurable space, let p, G M(X), let (V, +, •; t) be a 
locally convex space and let f : X —> V . we say that x £ V is the integral of f in 
X , denoted by 

(39) x:= f fdp 

J x 

if 

(40) = / * o fdp, 

Jx 

for all V £V*. 

It is known that if we have a compact metric space X, V is a locally convex 
space and / : X V is a continuous function such that cof{X) is compact then 
the integral of / in X exists and belongs to cof(X). We will also use the following 
well-known result, the barycentric formula: 

Proposition 3. [22] Let V be a locally convex space, let E C V be a complete, 
convex and bounded set, and p £ M 1 (E). Then there is a unique x £ E such that 




for all I £ V* . 

In the context of QIFS, we can take V = E = Mn- 
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6. Example: density matrices 



In this section we briefly review how the constructions of the previous section 
adjust to the case of density matrices. Define V :— Wn, V + := VHn (note that 
such space is a convex cone) , and let the partial order < on VH n be p < ip if and 
only if ip — p > 0, i.e., if ip — p is positive. Then 



is a regular state space [21] . Also, the set B of unity trace in V + is, of course, the 
space of density matrices, so B = A4n- 

Let Z C V* be a nonempty vector subspace of V* . The smallest topology in 
V such that every functional defined in Z is continuous on that topology, denoted 
by a(V,Z), turns V into a locally convex space. In particular, a(V, V*) is the 
weak topology in V. If (V, \\ ■ ||) is a normed space, then a{V* ,V) is called a 
weak* topology in V* (we identify V with a subspace of V**). We also have that 
(C, r) = (VHn,t), where r is the weak* topology (and which is equal to the 
Euclidean, see |H]) is a metrizable compact structure. In this case we have that 
B c = BHC = M N . 



Definition 4. A Markov operator for probability measures is an operator P : 
M 1 (X) -> M^X) such that 



forp 1 ,p 2 eM 1 (X), Ae (0,1). 

An example of such operator is the one which we have defined before and we 
denote it by V : M x (Mn) -> M X {M N ), 



We call it the Markov operator induced by the QIFS T = {Mn, Fi,Pi}i=i....,k- 
Define 

mb(X) :— {/ : X — > R : f is bounded, measurable} 
Then define U : m b (X) — > m b (X), 



{V,V+,e) = (n N ,Vn N ,tr), 



(41) 



p(Xm + (i - X)p 2 ) = ap^i + (i - X)PpL 2 




k 



(42) 



A- 



(43) 




Proposition 4. [21] Let f € rrib(X) and p € M 1 



(X) 7 then 




k 



(44) 



where (f,p) denotes the integral of f with respect to p. 



Definition 5. An operator Q : V + — > V + is submarkovian if 

(1) Q(x + y) = Q(x)+Q(y) 

(2) Q(ax) = aQ(x) 
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(3) \\Q(x)\\ < \\x\\, 

for all x, y G V + , a > 0. 

Every submarkovian operator Q : V + — > V + can be extended in a unique way 
to a positive linear contraction on V, see |21j . 

Definition 6. Let P : V + — > F + a Markov operator and let Pi : V + — > V^ + , i = 
1, . . . , fc 6e submarkovian operators such that P — J^. P^. We say £/ia£ (P, {Pi}i=i) 
is a Markov pair. 

From |21j . we know that there is a 1-1 correspondence between homogeneous 
IFS and Markov pairs. 



7. Some lemmas for IFS 
We want to understand the structure of A : TWjv — > -Mn, 

V iP V* 



(45) A(p) : y^p.P y^lnW.tM; ; ; 

i=l i=l \ i ) 

where 14, Wj are linear, J^. W*Wi = I. Such operator is associated in a natural way 
to an IFS which is not homogeneous. In this section we state a few useful properties 
which are relevant for our study. The following lemmas hold for any IFS, except 
for lemma El where the proof presented here is valid only for homogeneous IFS. 

Lemma 1. Let {X, Pj,Pi}i=i,...,fc be an IFS, \& a linear functional on X. Then 
U o = \J> o A, where U is given by ygp . 

Proof We have 

= £>(aO*(Fi(a;)) = *(^>(:r)P(x)) = *(A(z)) 



□ 



Corollary 1. Let T = (AT, P i ,p i ) i=lj ^ be an IFS and let p G X. Then A(p ) 
Po if and only if U(^(po)) = ty(p ), for all '5 linear functional. 

Proof Suppose that C(po) = pa- Then 

Utpipo)) = £>(p )*(Pi(po)) = *£>(po)Pi(po)) = *(A(po)) = *(po) 

Conversely, if W(*(p )) = *(Po), then 

*(A(p Q )) = W(*(p )) = *(po) 



□ 



Lemma 2. Le£ J 7 = {X, Fj,Pi}j=i,...,fc &e an IPS'. 

(1) Le£ p G A" such that Pj(po) — Po> * = 1, . . . , fc. T/ien Vop = 5 po . 

(2) Lei po G A" suc/i i/iai V<5 P0 = (5^,,, iften A(po) = Po- 
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Proof 1. We have 

VS P0 (B) = £ / PidS P0 = E / Pi(p)lB(Fi(p))dS p 

k k 

= ^2Pi(Po)lB(Fi(p )) = Epi(po)ls(po) = S P0 (B) 

i=l i=l 

2. Let \1/ be a linear functional. Then 



*(A(p )) = W(*(po)) = y = / *(p)dV<5 p 

*(p)d5po = *(po) 



□ 

Lemma 3. Lei {X, -Fi,Pi}i=i,...,fe &e a homogeneous IFS, A = ^Z-tPii^. 

(1) Lei &e £/ie barycenter of a probability measure v. Then A(p„) is the 
barycenter of Vv, where V is the associated Markov operator. 

(2) Let p be an invariant probability measure for V. Then the barycenter of p, 
denoted by p^, is a fixed point of A. 

Proof 1. We have, for ^> linear functional, 

*(A(p„)) = J ${A(p))dv = J Uo^dv = J VdVv 

2. By lemma ([1]), we have 

*(A(p M )) - U o 9{p M ) = Juo^dp = J VdVfi = J^dp = 
where the fact that U o $ is linear follows from the homogeneity of T . 

□ 

In order to prove uniqueness in item (2) above it would be necessary to assume 
hyperbolicity [3D]. It is known that without this hypothesis even in the classical 
case (for transformations for instance) it can happen the phenomena of phase tran- 
sition (two or more probabilities which are solutions) [23] [15] . The present setting 
contains the classical case and therefore in general there is no uniqueness. 



Example 3. Let k = N = 2, 



1 J 1 2 [ -3yJ 



3V2 
4 



Wi = (1/2)/, W 2 = (\/3/2)7. Then 



Hp) = J2p^(p) = J2tr(W iP W*) 



v iP v* 



tr(V lP V*) 
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induces an IFS and it is such that po = ^|0)(0| + ||1)(1| is a fixed point, with 
F\(po) — ^(po) = Pa- We can apply lemma [5] and conclude that 5 PQ is an invariant 
measure for the Markov operator V associated to the IFS determined by pi and Fi . 



The following lemma, a simple variation from results seen in |21) . specifies a 
condition we need in order to obtain a fixed point for A from a certain measure 
which is invariant for the Markov operator V. 

Lemma 4. Let {Mn, Fi,Pi}i=i t ...,k be an IFS which admits an attractive invariant 
measure p for V. Then rim„_ i . 00 A n (po) — p P , for every po G Mn, where p p is the 
barycenter of p. 

Proof Let p G Mn- Then 

*(A"( P0 )) = W B (*(po)) = J U n (V(p))d5 P0 = J V(p)dV n S P0 

so \I/(A n (po)) j x ^(p)dp = ^(pn), as n — > oo, for all linear functional. Hence, 
A™(p ) — > Pp. as n — > oo, for all p G Mn- 

□ 

In lemma 21 we have a general QIFS and an attractive invariant p, then p is 
the unique invariant measure, an easy consequence of attractivity [21] . In general, 
we will be interested in QIFS which has an attractive invariant measure. This will 
follow if we assume hyperbolicity. 

8. Integral formulae for the entropy of IFS 

Part of the results we present here in this section are variations of results pre- 
sented in |21) . Let (X, d) be a complete separable metric space. Let (V,V + ,e) 
be a complete state space, B = {x G V + : e(x) = 1} and T = {X,Fi,pi) i= i^,^k 
the homogeneous IFS induced by the Markov pair (A, {Ai}^ =1 ). Now define Ik := 
{!,..., fc}. Let n G N, I G 1%, i G h- Define F u := F. t o F L and 



(46) p Li {x) 



Pi(F L x)p L (x) if p L (x) ^ 
otherwise 



Proposition 5. Let n G N, / G rrib(X), x G X . Then 

{U n f)(x)=Y J P^)f(F i {x)) 

Proposition 6. Let x G B, n G N. Then 

A n (x)=J2p^) F ^)- 

Proposition 7. Let T be an IFS and let g : B — !> R. Then for n G N, 

(1) If g is concave (resp. convex, affine) thenU n g < goA n (resp. U n g > goA n , 
U n g = goA n ). 
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(2) If x is a fixed point for A then the sequence (U n g)(x)) n£ fi is decreasing 
(resp. increasing, constant) if g is concave (resp. convex, affine). 

Also suppose that T is homogeneous. Then 

(3) If g is concave (resp. convex, affine), then Ug is concave (resp. convex, 
affine). 



We recall some well-known definitions and results. Define ij : R + — > M. as 
, , . f — xloga; if x 4 

(47) ^ = {o if. = 

Then the Shannon-Boltzmann entropy function is h : X — ¥ R + , 

k 

(48) h(x):=J2v(Mx)) 

i=l 

Let n € N. Define the partial entropy H n : X —> M + as 

(49) H n (x) := *?(?*(*)), 

for n > 1 and Hq(x) := 0, x € X. Define, for x e X, 

(50) H(x) := limsup —H n (x), , %(x) := liminf —H n (x), 

rn-oo Tl n— >oo ft 

the upper and lower entropy on x. If such limits are equal, we call its common 
value the entropy on x, denoted by H(x). 

Denote by M V (X) the set of V-invariant probability measures on X. Let /i E 
M V (X). The partial entropy of the measure /i is defined by 

(51) H n (j*):= X>«P»>/*». 
for n > 1 and Ho(fi) '■= 0. 

Proposition 8. Let// <E M v (Jf). T7ien t/ie sequences (iiJ„(/i)) ne N and (Hn+iin) — 
Hn(p))neN are nonnegative, decreasing, and have the same limit. 

We denote the common limit of the sequences mentioned in the proposition above 
as and we call it the entropy of the measure [i, i.e., 

(52) H(fi) := lim -H n (fi) = lim (H n+1 (ji) - H n (fJ,)) 

n— >-oo 7i n— >oo 



The following result gives us an integral formula for entropy, and also a relation 
between the entropies defined before. We write S(fi) := M V (X) n Lim(V n //)„gN, 
where Lim( V™/i) n gN is the convex hull of the set of accumulation points of ( V™//) n gNj 
and Sj^(n) is the set iS(/t) associated to the Markov operator induced by the IFS 
J ' . For the definition of compact structure and (C, r)-continuity, see [2"T] . 
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Theorem 2. |21] (Integral formula for entropy of homogeneous IFS, compact case) . 
Let (C, r) be a metrizable compact structure (V,V + ,e) such that (A, {Aj}| =1 ) is 
(C, t)- continuous. Assume that po £ Bq ■= B f~l C is such that A(po) — Po- Then 

H(po) = U{v) = [ hdv 
Jx 

for each v £ Sjr c (6 po ), where Tc is the IFS T restricted to {Be, t). 

The analogous result for hyperbolic IFS is the following. 

Theorem 3. 21 Let T = (X, Fi,pi)i=i k be a hyperbolic IFS, x £ X, p, £ 
M 1 (X) an attractive invariant measure for J- . Then 

U{x) = lim (H n+1 (x) - H n (x)) 

n— >oo 

and 

H(x) = U(ji) = [ hdfi. 

Jx 

9. Some calculations on entropy 

Let U be a unitary matrix of order mn acting on H m <8 Hn- Its Schmidt decom- 
position is 

K 

U = Y^^V l A ®V l B i K = min{m 2 ,n 2 } 
»=l 

The operators V^ 4 and act on certain Hilbert spaces H m and H n , respectively. 
We also have that J2iLi 9i = 1- Let a — pa® p 1 ^ = Pa® In/n and define 

if 

Afcu) := tr B (UaU*) = ^q i V i A p A V i A * 

i=l 

Above, recall that the partial trace is 

tr B (\ai)(a 2 \ ® |6i)(6 2 |) := |ai)(a 2 |tr(|6 1 )(6 2 |) 

where |ai) and |a 2 ) are vectors on the state space of A and \b\} and |6 2 ) are vectors 
on the state space of B. The trace on the right side is the usual trace on B. A 
calculation shows that if p A = I m /m, then A(p A ) = p A and so A is such that 
A(J TO /m) = I m /m and A is trace preserving. 

Let T be the homogeneous IFS associated to the V A , that is, Pi(p) = tr(qiV A pV A *), 
F l {p) = (q t V A pV A *)/tr{q. l V A pV A *) and let p be a fixed point of A = 
Following [21], we have that p is the barycenter of V n S po , n £ N. By theorem [2 
we can calculate the entropy of such IFS. In this case we have 

(53) U(po)=H(v)= f hdv, 

J M N 

where v £ M V (X) n Lim(V ,l <5 Po )„ eN - 
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Let T = (Mn, Fi,pi)i = i,... y k be an IFS, A(p) = J2iPi-^i- Let U be the conjugate 
of V. By proposition El 

(U n h)(p) = J2 P'(p)KK(p)) 



N, 



and since h(p) = Ylj=i v(Pj(p))i we have, for i = (ii, . . . , i n ), and every po € M. 

(54) f hdV n S P0 = j U n hd8 po 

J Mn J Mn 

(55) = - / P'OO I>(*iO0) log Pj (F L (p))d5 P0 

JM N teJ n j- =1 



(56) = - P,(Po)Y,P^(Po))^m(F,(po)) 

(57) = - J] MPoK^po) ' ' •Pi„(^- 1 (^„_ 3 • • (^iiPo)))); 



(58) x ^p J (F I „(^„_ 1 (...(F I1 po))))logp J (F l , i (^ n _ 1 (...(^ 1 po)))) = {U n h){p ) 

3=1 

Suppose A(p ) = Pa- We have by proposition!?! since h is concave, that (U n h) n ^ 
is decreasing, W n /i < h o A™ and so 

(59) f hdV n S P0 < h(A n ( Po )) = h( Po ), 

J Mn 



for every n. 



10. An expression for a stationary entropy 



In this section we present a definition of entropy which captures a stationary 
behavior. Let H be a hermitian operator and Vi, i = 1, . . . , h linear operators. We 
can define the dynamics Fi : Mn — > Mn- 

VipV* 



(60) Fi(fi) 



tr{V iP V*) 



Let Wi, i = 1, . . . , k be linear and such that X)i=i W*Wi = I. This determines 
functions pi : Mn —> R, 

(61) Pi{p) := tr(W iP W*) 

Then we have X^iftG ) = L f° r ever y P- Therefore a family VF := {Wi}j=i,...,fc 
determines a QIFS J 7 ^/ = {Mn, Fi,pi}i = i....,k, with i^, pi given by ([60]) and (fUTj) . 
We introduce the following definition. 
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Definition 7. Let Tw be a QIFS such that there is a unique attractive invariant 
measure for the associated Markov operator V. Let pw be the bary center of such 
measure. Define the QIFS entropy; 
(62) 

HV{W) :=-X>K^/^)X>( ) logtr( 

I — 1 J — 1 

Remember that by lemma [4j we have that pw is a fixed point for 

k k 

(63) A(p)=Ar w (p):=J2pi(p)F i (p) = J2tr(W iP W*) 1 



tr(V iP V*) 

Lemma 5. hy{W) > 0, for every family Wi of linear operators satisfying W*Wi = 
L. 

Proof Note that, by definition, 

h v (W) = (Uh)(p w ) = [ hdV8 Pw 

J Mn 

and the function h (Shannon-Boltzmann entropy) is > 0. This proves the lemma. 
Another elementary proof is the following. Since pw is positive, we have that 
(W % p w W*v,v) = (p w W*v,W*v) > 0, v e Hn- So for {«j}j=i,...jv an orthonormal 
base for Hn, 

N 

tr{W iPw W*) = Y^WiPwWfmM) > 
i=i 

Analogously the expression above holds for the VipwV*, and therefore also for 
WjViPwVfW?, because 

{W 3 V lPw V*W*v,v) = (V iP wV*W;v,W*v) > 

To conclude that h v (W) > 0, we have to show that tr(WjVip w V*W*) < tr{V iPw V*). 
From W*Wi = I, we get 

k 

tr{W 3 V lPw V*W*) = tr{W*W 3 V lPw V*) < J2 tr ( W J W 3 V iPwV*) 

3=1 

k 

= tr (J2 WjWjVipwV?) = tr(V iPw V*) 

□ 

Remark For any fixed dynamics V, if we have that W^W m = I for some m 
then the remaining pi must be zero, because of the condition ^ W*Wi — I. In 
this case we have hy(W) = 0. We also have that hy(W) < logfc and for any given 
dynamics V, hviW) attains the maximum if we choose Wi — l/s/kl, for each i, 
where I denotes the identity operator. 
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Note that by the calculations made in section [9j we have hy{W) — Uh(pw), 
where Uh(p) = Y,iPi(p)H F i(p))- 

Lemma 6. Let T = (A4n, Fi,pi) be a QIFS, with Fj ; pi in the form i60\) and \61\) . 
Suppose there is po G Mn such that S Po is the unique V -invariant measure. Then 
A-f(pq) — po (Ajr is the operator associated to F) and 



U n hdS Po =U n h( Po ) = h(p ), 

for all n € N. Besides, U n h(po) — Uh(po) and so 

h v (W)=U n h(p ), 

for all neN. 

Proof The fact that A(po) = Po follows from lemma EJ item 2. Also, 
U n h( Pa ) = ju n M5 P0 = fhdV n S P0 = fhdS P0 = h(p ) 



and 

U n h{ Po ) = I U n hd5 P0 = I hdV n S P0 = I hdVS Po = I Uhd5 Pa = Uh{p ) 



□ 



Lemma 7. Let p be a V -invariant attractive measure. Then if p p is the barycenter 
of p we have, for any p, 

(64) lim U n h(p) = J Uhdp = J hdp < h(p p ) 

Proof The inequality follows from [3T], proposition 1.15. Also, by proposition 0] 
we have 

lim U n h{p) = lim / U n hdS p = lim / UhdV n - l 5 p = [ Uhdp, 

n— foo n— >oo J n— foo / ' / 

the last equality being true because of the weak convergence of (V n <5 p )„ e N- This 
proves the first equality in (|64|) . Since J Uhdp — J hdVp = J hdp, we obtain the 
second equality. 

□ 

Lemma 8. Let J 7 = (Aljy ,Fi,pi) be a QIFS, with Fi, in the form k6U\) and i61)) . 
Suppose that p is the unique point such that Ajr(p) = p. Suppose that Fi(p) = p, 
i = 1 , . . . , k . Then 

U n h{p) = h(p), 
n = 1,2, ... , and therefore hy(W) does not depend on n. 

Proof The proof follows by induction. Let n — 1. We have: 

uh(p) = J2Pi(pMFi(p)) - Hp) £>(p) - Hp) 

i i 

And note that U n h(p) = U(U n ~ 1 h)(p), which concludes the proof. 

□ 
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11. Entropy and Markov chains 

Let Vi, Wi be linear operators, i = 1, . . . , k, W*Wi = I. Suppose the Vi are 

fixed and that they determine a dynamics given by Fj : A4n — >• -Mn, i = 1, • • • , k. 
Define 

k 

(65) P := {(pi, . . . ,p fc ) : p; : .Mjv -> = 1, . • . , k, ^Pi{p) = 1, Vp G M»} 

i=l 

P' := P n {(px, . . . ,p k ) : 3Wi,i = 1, . . . ,k : Pi (p) = tr{W iP W*), 

(66) W l linear , ^ W^W; = /} 



(67) Af F := {M € M^jVIjv) : 3p e P' such that V p p = 
where V p : M 1 {M N ) -> M^Mjv), 

(68) V p (/i)(B) := E / ^ 

Note that a family W := {Wj}i = i j ... i fe determines a QIFS Tw, 

Let P = (Pij)i,j=i,...,N be a stochastic, irreducible matrix. Let p be the stationary 
vector of P. The entropy of P is defined as 

v 

(69) H(P) := - ]T PiPtj log PlJ 

We consider a few examples which will be useful later in this work. 
Example 4. (Homogeneous case, 4 matrices). Let TV = 2, fc = 4 and 



7 ' V 



/pw o y v >/pir 

Note that 

Poo + Pio 
Poi+Pn 

and so V*Vi = I if wc suppose that 



P := 



Poo Poi 
Pio Pll 



is column-stochastic. We have 



vipv; = ( P00 n Pl ° ) , v 2P v 2 * ( P0lPi 



J ' ' A \ 

w . i \ ) , vw = ( ° 

■ piopi J V PnPi 
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SO 

tr(VxpV*) = poopi, tr{V 2 pVi) = P01P4 
tr{V z pV^) =pi pi, tr{VipV£) =PnP4 



The fixed point of A(p) = J2i v iP v * is 







m , — 1-poo+Poi , 
\ 1— Poo+Poi / 

Let 7r = (tti, 7T2) such that Pit = it. We know that 

(70) *=( ^ , l ~ P ™ ) 

1 - Poo + Poi 1 - Poo + Poi 

Then the nonzero entries of pv are the entries of ir and so we associate the fixed 

point of P to the fixed point of a certain A in a natural way. Let us calculate 

hy(W). Note that A defined above is associated to a homogeneous IFS. Then 

Wi = Vi, i = 1, . . . , k and 

h v {W) = h v (V) 



(71) = - tr(v 3 V lPv V*V*) log (- 



tr(V lPv V* 



A simple calculation yields H(P) = hy(V), where H(P) is the entropy of P, given 
by ([69]) . This shows that the entropy of Markov chains is a particular case of the 
QIFS entropy. 







Example 5. (Nonhomogeneous case, 4 matrices). Let N = 2, k = 4 and 
Vi= ( JpTo \ ^ _ / ^TT 



/ ' V 



PlO J \ y'pTT 



Note that 



/' V 

n- / °\ ( 



Poo+Pio \ v^tt^ti/ ^ 9oo + 9io 



poi+Pn / \ 901+911 

and so J2i Vi*Vi = J+ W*Wi — I if we suppose that 



P:=( 


' Poo 


Poi A 


,<?:=( 


' 9oo 


9oi A 




V PlO 


Pn / 




v 9io 


9n / 
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are column-stochastic. Then 

tr(VipV*) = P00P1, tr{V 2 pV£) = p m p± 
tr(V 3 pV 3 *) = pwPi, tr(VipVl) = p n p 4 
tr(WipW*) = qoopi, triyVipW^) = q ip4 
tr{W 3P W;) = q w pu tr{W iP Wl) = qup 4 
We want the fixed point of A(p) — ^Z i Pi{p)F i {p). This leads us to 

goo ( P00P1 \ %i ( pmPi \ qw_ ( \ On / 
Poo V / poi V / p 10 V PioPi / Pn V piip 4 



P 



Note that the Pij cancel and so we obtain a calculation which is the same as the 
one obtained in the previous example. Hence 

\ 



nixr — l-9oo+9oi | 
P w — I n 1-900 I ' 

\ l-9oo+9oi / 

and its nonzero entries are the entries of the fixed point for the stochastic matrix 
Q. Calculating h v (W) gives 

I 1 J 1 

(72) 

-(<?oolog<7oo + <?iologgio) — — (qoi log 901 + 9n log 911) = H(Q) 



qoi +9io " 9oi + 9io 

So we have obtained a calculation which is analogous to the one for the homogeneous 
case. This result generalizes what we have seen in the previous example. 







Example 6. (Homogeneous case, 2 matrices). Let N = 2, k = 2 and 
Vl = ( ^ ° V V 2 - ( ° ^ 



/Pw / V y^pIT 

Note that, just as in the previous examples 



i ^ 

and so J2i Vi^i = ^ if wc suppose 

P := 

is column-stochastic. The fixed point for A is 

Pv = 



Poo + Pio 
P01+P11 



Poo Poi 
Pio P11 



Poi P00P10P01 1 P01P11P10 

P01+P10 P01+P10 P01+P10 
P00P10P01 1 poipupio pio 
poi+pio poi+pio poi+pio 
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The entries of the main diagonal of pv correspond to the entries of the fixed point 
of P. The entries of the secondary diagonal are a linear combination of the ones in 
the main diagonal. Then for the Vi chosen we have 

MVjV iPv v?v;y 



(73) h v (W) = h v (V) = -Y J tr(y j V i pvV'V*') log (- 



• ■J 



tr{V lPv V*) 



H(P) 



by an identical calculation made for the equation (|72|) from the previous example. 
In other words, the fact that the fixed point of A is not diagonal does not change 
the calculations for the entropy. 



Example 7. (Nonhomogeneous case, 2 matrices). Let N — 2, k = 2, 



Vx 



Wi = 



/Poo 

/Pw 

/9oo 

/9io 



V 2 



Wo = 



/Poi 
/PIT 

Von 



As in the other examples, J2i — Yli W*Wi = I if we suppose 



P:=| 


' Poo 




Q:= ( 


' 9oo 


9oi \ 




V Pw 


pu J 




v 9io 


9n / 



is column-stochastic. From 

tr{V lP V?) = Pu tr(V 2 pV 2 *) = p 4 
tr(W lP W;) = p u tr(W 2P WZ) = Pi 
tr{WxVi P V;Wl) = poopi, tr(W 2 V lP V*WZ) = p wP i 
tr{W x V 2P V^Wi) = poiPi, tr(W 2 V 2 pV 2 *W 2 *) = p n p 4 
and a simple calculation, we get hy(W) — H(P). 



Lemma 9. Let Vij be matrices of order n, 



for i,j = 1, 



,n. Let 



where P = (pij)ij=i,. 

Proof Note that 
(74) 



so 



A|(p) 



/Pij\i){j\ 



1,3 



Then for all n, A n P (p) = A P n(p). 
VkiVi 3 = ^/pki^/pljSu\k)(j\ 



= E E I fc > tflPlj) W = E Phi \k) (J\P\J) (k\ = A P2 (p) 

k,j i k,j 

The general case follows by iterating the above calculation. 



A THERMODYNAMIC FORMALISM FOR DENSITY MATRICES IN Q. I. 



25 



□ 

Corollary 2. Under the lemma hypothesis, we have limn-^ A 7 p(p) = A 7T {p), where 
7r = lim n _ i . 00 P n is the stochastic matrix which has all columns equal to the station- 
ary vector for P. 

12. CAPACITY-COST FUNCTION AND PRESSURE 

Recall that every trace preserving, completely positive (CP) mapping can be 
written in the Stinespring-Kraus form, 

fe fe 

A(p) = J2 v ip v *> E^ = / < 

»=1 i=l 

for Vi linear operators. These mappings are also called quantum channels. This is 
one of the main motivations for considering the class of operators (a generalization 
of the above ones) described in the present paper. These are natural objets in the 
analysis of certain problems in quantum computing. 

Definition 8. The Holevo capacity for sending classic information via a quan- 
tum channel A is defined as 

n n 

(75) C A := max s( J^PiHPi)) ~ I>s( A (pi)) 

PiGM N i=1 i=1 

where S(p) — —tr(p\og p) is the von Neumann entropy. The maximum is, therefore, 
over all choices of pi, i = 1, . . . ,n and density operators pi, for some neff. The 
Holevo capacity establishes an upper bound on the amount of information that a 
quantum system contains [17] . 

Definition 9. Let A be a quantum channel. Define the minimum output en- 
tropy as 

H mm (A) :=mmS(A(|V>}(V>l)) 

\4>) 

Additivity conjecture We have that 

Ca 1 ®a 2 = Cai + Ca 2 
Minimum output entropy conjecture For any channels Ai and A2, 
iT nm (A! ® A 2 ) = i? mm (Ai) + H min {A 2 ) 

In [19], is it shown that the additivity conjecture is equivalent to the minimum 
output entropy conjecture, and in |10) a counterexample is obtained for this last 
conjecture. 

Remark Concerning QIFS, our interest in capacity is motivated by the following 
observation. Considering expression (I75|) . note that the term 

n 

(76) 5>S(A( ft )) 

i=l 
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is a convex combination of von Neumann entropies, in the same way as the QIFS 
entropy. So we see that given a QIFS, we can consider capacity functions, and the 
QIFS entropy arises in a natural way. For an example, we perform the following 
calculation. If A^ are the eigenvalues of p then we can write 

(77) S{ P ) = -]T A, log \ 

i 

Then write the QIFS entropy as 

k 

(78) h v (W) = -J2 tr{W i pwW*)ai j (pw)loga ij (pw) 

*>i=i 

where 

(79) aijip) := tr(V iP V*) 

We see that for p w e Xjv and i fixed, we have ^ ciij(pw) — 1- Define for each i 
the density operator 

(80) A :=$^Oij(pw)|j)0'l 

3 

Then by ([77]), 

(81) S(pi) = ~y^ j ajj (p w ) log a-ij (pw) 

3 

By (|78|) . we can write 

(82) h v {W)( Pw ) =Y / tr(W i pwWl)S( Pi ) 

i=l 

A Positive Operator- Valued Measurement (POVM) is described by a set of 
positive operators Pi (POVM elements) such that J^i Pi = I- If the measurement 
is performed on a system described by the state vector |^), then the probability of 
obtaining i as the outcome is given by 

(83) Pi = 

Note that a QIFS F induced by linear V* and W l7 contains a POVM by taking W*W l 
as POVM elements. If X is a random variable that takes values p±, . . .pk then the 
Shannon entropy is H(X) = — Y2iPi ^°SPi an d the joint entropy of variables X and 
Y is 

(84) H(X,Y) := -J2p&v)logp(x,y) 

where p(x, y) is the probability that X = x and X = y. The mutual information 
I(X : Y) is defined by I(X : Y) := H(X)+H(Y)- H(X,Y). Then, considering the 
QIFS entropy we can state the Holevo bound in the following way: first consider a 
QIFS J- such that there is a unique attractive measure which is invariant for the 
Markov operator V associated to T . Let pw be the barycenter of such measure. 



A THERMODYNAMIC FORMALISM FOR DENSITY MATRICES IN Q. I. 



27 



Theorem 4. (Holevo bound for QIFS) Suppose J- is induced by linear operators 
Vi and Wi with J2i W*Wi = I and for each i = 1, . . . , k write pi — tr(WipwW*) 
and pi = J2j a ij(Pw)\j)(j\, where ay is given by J7ff| ). Suppose Alice prepares a 
mixed state px chosen from the ensemble {p\, ...,pk} with probabilities {p\, . . . ,Pk} 
(that is, we assume px is a state determined by a random variable X such that it 
assumes the value pi with probability pi). Suppose Bob performs a POVM measure- 
ment on that state with POVM elements {-Pi}i=i,..., m and measurement outcome 
described by a random variable Y . Then, by writing p = ^iPiPi, we have 

(85) I{X:Y)<S{p)-h v {W){ Pw ) =:(,{£) 

The number £(£) is the Holevo information of the ensemble given by £ = 
{PilPi}i=i,...,k- We see that ([55)) holds by applying the Holevo bound for the von 
Neumann entropy (see [17]) together with ([80[) and ([82[) . 



We are also interested in a different class of problems which concern maximiza- 
tion (and not minimization) of entropy plus a given potential (a cost) [5].|11].[I2]. 

Definition 10. Let Mp be the set of invariant measures defined in the section \H\ 
and let H be a hermitian operator. For p £ A4p let p^ be its barycenter. Define 
the capacity-cost function C : M + — ► M + as 

(86) C(a) := max {hwyifiu) ■ tr(Hpn) < a} 

The following analysis is inspired in [8j. There is a relation between the cost- 
capacity function and the variational problem for pressure. In fact, let F : R + —> R + 
be the function given by 

(87) F(X) := sup {h W y{p^)-\tr{Hp^)} 

fj.eM F 

We have the following fact. There is a unique probability measure vo £ Mf such 
that 

F(X) = hw,v(p» ) - \tr{Hp Va ) 
Also, we have the following lemma: 

Lemma 10. Let A < 0, and a = tr(Hp Vo ). Then 

(88) C(a) = h wy ( Puo ) 

Proof Let v € A4f, v ^ vq, with tr(Hp v ) < a = tr(Hp Va ). Then 
hw,v(Pv) ~ \tr(Hp v ) < hw,v{pu ) ~ Xtr(Hp Va ) 

so 

hw,v{pu) < h W y(p Uo ) 

Hence 

hw,v{Pv a ) = sup {hw,v(Pn) ■ tr(Hp^) < a} = C(a) 

□ 
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13. Analysis of the pressure problem 
Let Vi, Wi, Hi be linear operators, i — 1, . . . , k, with J2i W*Wi — I and let 

fc 

(89) Hp:=J2 H iP H t 

i=l 

a hermitian operator. We are interested in obtaining a version of the variational 
principle of pressure for our context. We will see that the pressure will be a max- 
imum whenever we have a certain relation between the potential H and the prob- 
ability distribution considered (represented here by the Wi). We begin by fixing a 
dynamics, given by the Vi. From the reasoning described below, it will be natural to 
consider as definition of pressure the maximization among the possible stationary 
Wi of the expression 

fc 

h v (W) + lo S (tr{H 3 ppH*)tr{V 3P pV*))tr{W ]Pw W*) 

3=1 

where pp is the eigenstate of a certain Ruelle operator, described below. We begin 
our analysis by using the following elementary lemma. 

Lemma 11. |18j If rx, . . . , and th, . . . , are two probability distributions over 
1, . . . , k, such that r 3 > 0, j = 1, . . . , fc, then 

fc fc 

(90) - J2 a 3 lo § U + q J lo § r 3 - 

i=i i=i 

and equality holds if and only if rj = qj, j = 1, . . . , fc. 



The potential given by (|89)) . together with the V i7 induces an operator given by 

k- 

(91) C„{p) -Y^triHipHDVpV* 

i=l 

By proposition [5] we know that such operator admits an eigenvalue /3 with its 
associated eigenstate pp. Then Ch(pp) = fipp implies 

k- 

(92) tr{Hip p Hi)ViPpV: = ppp 

i=l 

In coordinates, (1921) can be written as 

fc 

(93) ^trlHippHfttVippVnfl = P(pp)ji 

i=l 

Remark Comparing the above calculation with the problem of finding an eigen- 
value A of a matrix A — (ay), we have that equation (|92p can be seen as the 
analogous of the expression 

(94) IE A = XI 
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Above, the matrix A plays the role of a potential, E A denotes the matrix with 
entries e aij and lj denotes the j-th coordinate of the left eigenvector I associated 
to the eigenvalue A. In coordinates, 



From this point we can perform two calculations. First, considering (|92|) we will 
take the trace of such equation in order to obtain a scalar equation. In spite of 
the fact that taking the trace makes us lose part of the information given by the 
eigenvector equation, we are still able to obtain a version of what we will call the 
basic inequality, which can be seen as a QIFS version of the variational principle 
of pressure. However, there is an algebraic drawback to this approach, namely, 
that we will not be able to recover the classic variational problem as a particular 
case of such inequality (such disadvantage is a consequence of taking the trace, 
clearly). The second calculation begins at equation (|93|) . the coordinate equations 
associated to the matrix equation for the eigenvectors. In this case we also obtain 
a basic inequality, but then we will have the classic variational problem of pressure 
as a particular case. 

An important question which is of our interest, regarding both calculations men- 
tioned above, is to ask whether it is possible for a given system to attain its max- 
imum pressure. It is not clear that given any dynamics, we can obtain a measure 
reaching such a maximum. With respect to our context, we will remark a natural 
condition on the dynamics which allows us to determine expressions for the measure 
which maximizes the pressure. Now we perform the calculations mentioned above. 

Based on ([92]) . define 



(95) 








(96) 



tj = ~tr(H jP0 H*)tr(V jP pV;) 



So we have J2j r j = 1- Let 



(97) 




where, as before, pw is the fixed point associated to the operator Ajr, 



w 



k 



(98) 




i=i 



induced by the QIFS {Mn , Fi,Pi)i=i, 



Hp) 



v iP v* 



tr{V lP V*) 



and 



Pi(p) 



tr(W iP W*) 
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Note that we have 

k 

= wm^vn% w ! w ' v,pwV:} = 1 

Then we can apply lemma [TT] for rj, q l j, j = 1, . . . k, with i fixed, to obtain 



^ V tr(V,owV*) ) l0gtr \ 



tr{V lPw V*) ) b V tr{V lPw V* 



"(VSpwrV,*) V °V/3 

and equality holds if and only if for all 
(100) ~tr(H jP pH*)tr(V jPf) V*) 



Then 



v iPw v* w* \ / Wj^pw 



ir-(y iPlv K*) ) b \ tr{V tPw V*) 



tr{V lPw V % 



- ^ V tr(Vi ou/V*) J &H 



3 

which is equivalent to 

-W 3 V lPw V*W* x ^ , fWjVjpwVCW; 



( 101 ) + E fr^vr) los (M*^ <i°g/? 

Multiplying by tr(Wi P wW*) and summing over the i index, we have 

tr(W iPw W*) 



(102) <E^WWW7) log/3 = log/3 
and equality holds if and only if for all i, j, 

(103) ^(WW) - 1^4%^ 
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Let us rewrite inequality (|102|) . First we use the fact that pw is a fixed point of 

k 



(104) ^tr{W iPw W*) 



tr{V lPw V* 



Pw 



i=l 

Now we compose both sides of the equality above with the operator 

k 

(105) ^ log (triHjppHfltriVjppV^WfWj 
i=i 

and then we obtain 

k k 

tr{W iPw W*) J2 l°g (tr(H 3 p p H*)tr(V ]Pp V*)) W*W 3 

k 

(106) = Pw (tr(H jPf} H*)tr(V jP pV*)) W*W 3 - 

i=i 

Reordering terms we get 

± log (hh^mvjppv;)) ± t7 ^^:h pwv;w;w 3 

i — l \ r i 



3 = 1 



(107) = p w ]T log (tr(H jP pH*)tr(V jP0 V;)) W*W 3 

i=i 

Taking the trace on both sides we get 

k k 

£ log (tr(H jPp H*)tr{V jPfs V*j) £ tr( Wj V iPw V*W*) 



3 = 1 



(108) = E lQ g (tr{H 3 PpH*)tr{V jP pV*))tr{ Pw W*W 3 ) 

3 = 1 

Note that the left hand side of (|108[) is one of the sums appearing in (|102j) . Therefore 
replacing f| 108[) into (|102l) gives our main result. 

Theorem 5. Let Tw be a QIFS such that there is a unique attractive invariant 
measure for the associated Markov operator V. Let pw be the bary center of such 
measure and let pp be an eigenstate of Ch{p) with eigenvalue f3. Then 

k 

(109) h v (W) + ^ log (tr{H jP pH*)tr{V jPp V*))tr{W jPw W}) < log/3 

3 = 1 

and equality holds if and only if for all i,j, 



1 tr(WiVip w Vi*W?) 
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In section [16] we make some considerations about certain cases in which we can 
reach an equality in (|109[) . 



For the calculations regarding expression (|93l) , define 

(HI) r jlm = -tr(H j p p H j ) J 

P {Pp)lm 

Then we have Y]j rji m — 1 . Let 

(WjV,p w v-w; 



(112) *-n t< iwin 

A calculation similar to the one we have made for (|109[) gives us 

k 

h v (W)+Y, tr(Wj Pw W* ) log tr(Hj Pfj H* ) 



3=1 



(113) + tr(W jPw W*) log ( ™^ to ) < bg p 

and equality holds if and only if for all I, m, 

,„ A s 1 rrr tr(WjVjpw V? W* ) 
114 -tr{H lP pH J = — 







14. Revisiting the eigenvalue problem 
Consider the operator 

k 

(115) C H (p) = tr{HipH*)V iP V? 

i=l 

induced by a fixed dynamics Vi i = 1, . . . , k, Vi linear, and by Hp := ^ HipH* , 
i?j linear. The eigenvalues equation for Ch written in coordinates gives us the 
following system, for k = 2: 

tr(H lP pH*)(vl lP u + 2i>iivi2Pi2 + vf 2 p 22 ) 

(116) + tr(H 2 p H 2 )(wl 1 pu + 2wnwi 2 pi2 + w\ 2 p 22 ) = P>n 

tr(H 1P pHi)(v 2 iV 11Pll + (v 21 vi3 + V22Vn)p 12 + ^22^12^22) 

(117) + tr(H 2p/ 3H 2 )(w 21 Wi lP u + (w 2 iWi 2 + w 22 wu)p 12 + W22W12P22) = PP12 



(118) 



tr(Hi P ^H*)(vl lP n + 2V21V22P12 + v\ 2 p 22 ) 

+ tr(i? 2 P/3ff 2 *)( w 2lPll + 2w 2 lW22Pl2 + ^22^22) = P>22 
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And we can also write, for i = 1,2, 
(119) 

tr(H iP H*) = ((h^) 2 + (h\ 2 ) 2 )pu + 2(h\ x h\ 2 + h\ 2 h\ 2 ) Pl2 + ((h\ 2 ) 2 + (h\ 2 ) 2 )p 22 







Fix H U H 2 , let Vi, V 2 be defined by 

(120) Vl =( ^ 7), V 2 =( ° ° 

V J \ w 2 i w 22 

then we have, by (|116[) - (|118[) that p\ 2 — and 

(121) triHippHD^pn + v\ 2 p 22 ) = (3p n 

(122) tr(H 2 ppH*)(w 2 lPll + w\ 2 p 22 ) = pp 22 
that is, 

(123) [((h^) 2 + (h\ 2 ) 2 ) Pll + ((h\ 2 ) 2 + (h\ 2 ) 2 )p 22 ](v 2 llPll + v\ 2 p 22 ) = p Pll 

(124) [((h 2 ,) 2 + (h\ 2 f) Pll + ((h 2 2 f + (h 2 22 ) 2 )p 22 }(w 2 2lPll + w 2 22 p 22 ) = (3p 22 
Also, suppose that 

(125) V\x = v V 2 = w 2 i = w 22 = 1 
Then we get 

(126) ((h^) 2 + (hl 2 ) 2 ) Pll + ((hl 2 ) 2 + (h\ 2 ) 2 ) P22 = fj Pll 

(127) ((h 2 n ) 2 + (h 2 2 ) 2 ) Pll + ((h 2 2 ) 2 + (h 2 22 ) 2 ) P22 = (3p 22 



Let A — (cty) be a matrix with positive entries and consider the problem of 
finding its eigenvalues and eigenvectors. Then from 

(128) fflnui + a 12 v 2 = ftvi 

(129) a 2 \Vi + a 22 v 2 = f3v 2 

we see that the systems (|126|) - (|127p and (]128|) - (|129[) are the same if we choose 

(130) aii = (^ii) 2 + (^i2) 2 , aia = (^ 2 ) 2 + (h\ 2 ) 2 

(131) an = (h 2 n ) 2 + (h 2 12 ) 2 , a 22 = (h\ 2 f + (h\ 2 ) 2 

We conclude that Perron's classic eigenvalue problem is a particular case of the 
problem associated to Cr acting on matrices. In fact, if we fix 

(132) Vt=(l I), V 2 '° ° 



jo/' v 1 1 

and given A a matrix with positive entries, choose 
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Then the operator Lh has a diagonal eigenstate 



(134) 



PP = 



Pn 

P22 



associated to the eigenvalue /3, and we have that, defining v = (pn,p22)> we get 
Av = j3v. 



1 1 




Example 8. Let 

(135) Vi = 

Then Av = f3v leads us to 
(136) 



, \ A 1 4 



vi + Av2 — (3vi 



(137) 

The eigenvalues are 
with eigenvectors 



3wi + -v 2 = (3v 2 



4 4 



( _L ± J-VT93, 1) 



l±£ + ^/l93 v 12 12 



193, v= 7^ ( 1 ^ + TC\/I93 i 1) such that Av = fiv. 

J-TTTT To V J- Jo 1Z 1Z 



Then we have ^ = f + i- 777 ^ 1 ' 1 ' 1 

Let 
(138) 



Hi 



/all 
^/ali" 



1 
2 



#2 



/02T 
\/a22 



V3 

71 



Then solving Ch(p) = ftp gives us pi2 = and 
(139) pn + 4p 22 = Ppn 



(140) 



3P11 + = /fy>22 



which is the same system as (|136[) - (|137|) . So (3 = | + jVT93 and the corresponding 
eigenstate, since p\ 2 = 0, is 



'193 



(141) 







!+T5 + T2V193 
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15. Some classic inequality calculations 

A natural question is to ask whether the maximum among normalized Wi, i = 
l,...,k, for the pressure problem associated to a given potential is realized as the 
logarithm of the main eigenvalue of a certain Ruelle operator associated to the 
potential Hi, i — 1, . . . , k. This problem will be considered in this section and also 
in the next one. 

We begin by recalling a classic inequality. Consider 

k k 

(142) - Y Qj log Qi + X 9i lo S r o < 

3=1 3=1 

given by lemma 111! Let A be a matrix. If v denotes the left eigenvector of matrix 
E A (such that each entry is e oy ), then vE A — f3v can be written as 

(143) £>e°« =/3^-, Vj 

i 

Define 
(144) 



So J2i r ij = 1- Let lij > such that q^ = 1. By (|142p . we have 

(145) - J2 Qij log Qij + Y qi i log ~q^' - 

i=l i=l " i 

That is, 

k k k 

(146) ~Y Ho 1°S 3w + 51 a ^ + X qi i ( log Vi ~ log U J ) - log $ 

i—1 i—1 i—1 

Let Q be a matrix with entries q^, let n — (jri, . . . , 7Tfc) be the stationary vector 
associated to Q. Since J^iQij = 1, Q is column-stochastic so we write Qn = it. 
Multiplying the above inequality by 7Tj and summing the j index, we get 

(147) y ~< X : " g ' / " ■ X X ,// '"" • X t ' X qvQ°z r < l °g v j) ^ i°g^ 

j i j i j i 

In coordinates, Qix — 7r is ■ <7ij7Tj = 7^, for all i. Then 



X~ X *i log ''■> ' X X 



Qij Qij 



(148) + X ^ X *j log Vl ~ X ^ X ^ log W J - log 13 

j i j i 

These calculations are well-known and gives us the following inequality: 

(149) ~ X ^ X 9*3 log Qv + X n J X 

Qij O'ij <log/3 

j i j i 

Definition 11. We call inequality the classic inequality associated to the 

matrix A with positive entries, and stochastic matrix Q. 
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Definition 12. For fixed k, and I, m — 1, . . . , k we call the inequality 

k 

h v (W) + J2 tr{W oPw W*)\ogtr{H ]P pH*) 

3=1 

(150) + Y, tr{W 3 p w W*) log ( {Vj ^ V J )lm ) < log ft 

the basic inequality associated to the potential Hp = J^. HipH* and to the QIFS 
determined by Vi, Wi, i = 1, . . . , k. Equality holds if for all i, j,l,m, 



(151) -tr{H ]Pp H 



1 .„ ->{W7)lm _ trQVjViPwV^W;) 
ptr^ppHj) tr(V lPw V*) 





As before, pp is an eigenstate of £h(p) and pw is the barycenter of the unique 
attractive, invariant measure for the Markov operator V associated to the QIFS 
Fw- Given the classic inequality (|149l) we want to compare it to the basic inequality 
(|150[) . More precisely, we would like to obtain operators Vi that satisfy the following: 
given a matrix A with positive entries and a stochastic matrix Q, there are Hi and 
Wi such that inequality (|150[) becomes inequality (|149[) . We have the following 
proposition. 

Proposition 9. Define 

Let A = (ay) be a matrix with positive entries and Q — (qij) a two-dimensional 
column- stochastic matrix. Define 

< 153 » ""= )■*— ( 



>>'>-« = ! .Lr Lt). ",2 = 



/g a 21 */g&21 / \ A /pQ22 \/ p a 22 



and also 

(155) Wi 



( \/<7iT 


°) 


, w 2 = [ 


f 




v 







v o 


o J 



(156) 



W 3 = 



^ o • w < = ( n a- 

TTien the basic inequality associated to Wi,Vi,Hi, i = 1, ...,4, Z = m 
I = m = 2, is equivalent to the classic inequality associated to A and Q. 



1 or 



We use the following lemma. 
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Lemma 12. For Vi given by 



(157) 
(158) 



Vi = 



V,= 



fvTi 





/iii 



Vo = 



/Vi2 






where vij > 0, we have that the associated QIFS is such that pw cmd t 
density operators, for any choice of Wi and Hi, i = 1, . . . , 4. 

Proof of Lemma 1121 We have that pw is a fixed point of 

v iP v* 



are diagonal 



A(p)=J2HW iP W* 



tr{V lP V* 



Writing 



Pn Pl2 
Pl2 P22 



we have that A(p) = p leads us to 

tr{WipW*) ( vnpn 




tr(W 3P Wi 



tr(W2pW^) ( v 12 p 2 2 




tr(V 2 pV*) 







tr{V 3P V 3 * 





V2ipn 



tr{V iP Vi) 








H 


' Pll 


P12 \ 




U22P22 j 




\ P12 


P22 J 



Then pi2 — and so pw is diagonal. In a similar way we prove pp is diagonal. 



□ 



Proof of Proposition [9] Let Vi, Wi, i = 1, ... ,4 and Hij, i,j = 1, 2 as in the 

statement of the proposition. A simple calculation shows that 

(159) triHijppH^) = e^ 

(since pp is diagonal, by lemma [T2")) . By example O the choice of Vi and Wi we 
made is such that the entropy hv(W) reduces to the Markov chain entropy. Then 
a calculation yields 

(ViPeXQu _ {P )u 



(160) 
In a similar way, 
(161) 



{Pf})n (pp)n 

{ViPpV*) 2 2 _ (Pfi)22 



1 



(Pp)22 (p • 13)22 

Then from the basic inequality with I = m = 1 or 



2 we get 



(162) h Y (W) + tr(W jPw W*) ^nr ^Jr^ logtr(Hi P pH?) < log 



Finally, since tr(HijppH*j) 
(fT49l) . 



and Q77 = 77, we conclude that (|162p becomes 

□ 
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Example 9. Let 

Then 

If we suppose the Vi are the same as from proposition[9l we have that pp is diagonal, 
so 

tr(HippH*) = 4, tr{H 2Pfi H* 2 ) = 1, tr(H 3 p H* ) = 2, tr(H iPfj H*) = 1 
Then Cu(p) = /3/0 leads us to 

4/011 + P22 = /3pil 
2/011 + P22 = /3/022 

A simple calculation gives 



with cigcnstatc 



o 



7 + Vrf 1 1 





We want to calculate the Wi which maximize the basic inequality (|150l) . Recall 
that from proposition |9l the choice of Vi we made is such that 

So 

fe 

(163) /iy(W) + £ fr(Wj/>wrW7) \ogtr{H ]Pp H*) < log/3 

3=1 

and equality holds if and only if, for all i, j, I, m, 



(164) -tr{H ]Pp H 



Choose, for instance, I = m = 1. Then condition (|1 64[) becomes 

(165) ^(W f )= Mvwvy) 

To simplify calculations, write Wi = W*Wi and Wi = {w\j). Then we get 

(166) tr y =<=4, < = 1 4 

So we conclude 
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That is, 

(108) w, - w 2 - -L/, . . -L/ 

Note that 

To solve that, we renormalize the potential. Define 



(169) ffj := V^Ht, a := 



4 + V2 



Then a calculation shows that Cf I {p) — ftp gives us the same eigenstate as before, 
that is pp = pp. But note that the associated eigenvalue becomes $ — a(3. Now, 
note that it is possible to renormalize the Wi in such a way that we obtain Wi with 
^2iW*Wi = I, and that these maximize the basic inequality for the Hi initially 
fixed. In fact, given the renormalized Hi, define 

(170) Wi = y/aWi, * = 1,...,4 
Note that £\ W*W = I. Also we obtain 

k 

(171) h v (W) + J2tr(W jP ^W*)logtr(^H jP ^H*) < loga/3 
which is equivalent to 

k 

(172) h v (W) +J2 tr (WjPwW*) \og{atr{H jPf) H*)) <\oga + log/3 

3=1 

That is 

fe 

^(W) + ^tr(^p#W*)loga 

3=1 



k 

(173) +^tr(^p^W*)logir(^p^;) < loga + log/3, 

and if we cancel log a on both sides, we get the same inequality as for the non- 
renormalized Hi. As we have seen before, such Wi gives us equality. Hence 

k 

(174) h v (W) +J2tr(W j PyyW*)logtr(H j p^H*) = log/3 

3=1 
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16. Remarks on the problem of pressure and quantum mechanics 

One of the questions we are interested in is to understand how to formulate a 
variational principle for pressure in the context of quantum information theory. An 
appropriate combination of such theories could have as a starting point a relation 
between the inequality for positive numbers 

- q% log * + Y * lo S Pi - °' 

i i 

(lemma II H seen in certain proofs of the variational principle of pressure) , and the 
QIFS entropy. We have carried out such a plan and then we have obtained the 
basic inequality, which can be written as 

fc 

(175) hv(W) + Y {tr{H ]Pl3 H*)tr{V ]P pV;)yr{W 3 p w W*) < log/3 

3=1 

where equality holds if and only if for all i,j, 



(176) fr{H mHj )tr{V 3 p,V 3 ) = fr( } 



As we have discussed before, it is not clear that given any dynamics, we can obtain 
a measure such that we can reach the maximum value log /3. Considering particular 
cases we can suppose, for instance, that the Vi are unitary. In this way we combine 
in a natural way a problem of classic thermodynamics, with an evolution which has 
a quantum character. In this particular setting, we have for each i that ViV* — 
V*Vi = I and then the basic inequality becomes 

fc 

(177) h v (W) + Y tr(Wjp w W*)logtr(H jP pH*) < log/3 

3=1 

and equality holds if and only if for all i, j, 

(178) ^tr(H ]Pp H*) = tr{W 3 V lPw V*W;) 
We have the following: 

Lemma 13. Given a QIFS with a unitary dynamics (i.e., Vi is unitary for each 
i), there are Wi which maximize Ijl75\ l, i.e., such that 

fc 

(179) h v (W) + Y J HW J Py V W;)\ogtr{H p p H*) = log/3 

3 = 1 

Proof Define, for each j, 



(180) W j := yjjtriHjPpHfll 

where / is the identity. The equality condition (I178P is satisfied by such Wj , so the 
lemma follows. 

□ 
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Remark The above lemma also holds for the basic inequality in coordinates, 
given by (|150|) . Also, it is immediate to obtain a similar version of the above 
lemma for any QIFS such that the Vi are multiples of the identity, and also for 
QIFS such that pw fixes each branch of the QIFS, that is, satisfying, for each i, 

ViPwV* 
tr(V iPw V*) PW 



17. Concluding remarks 

Considering the QIFS setting, we defined a concept of entropy and a Ruelle 
operator in such a way that we are able to get some analogous results to the classical 
Thermodynamic Formalism. Such Ruelle operator admits a positive eigenvalue, 
which gives us an upper bound for the pressure (entropy plus a potential) associated 
to the QIFS. Our configuration space is the set of density matrices. We did not 
consider the usual space of symbols or a shift operator, as it is assumed in the 
Ruelle-Perron-Frobenius theory. We have replaced the dynamics given by the shift 
with the one given by the inverse branches of the iterated functions (which are 
defined by a set of operators). 

The references [16] and [21] are of fundamental importance in our investigation. 

A starting point for further investigation could be to study more properties of the 
QIFS entropy, such as convexity and subadditivity. Also, a natural question is to 
ask whether it is possible to consider a QIFS acting in an infinite tensor product of 
finite Hilbert spaces which would be the analogous of considering the full Bernoulli 
space. 

In a forthcoming paper we are going to consider relative entropies and quantum 
conditional expectations. 
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